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Dark soliton formation in mode-locked lasers is investigated by means of a power-energy saturation model 
which incorporates gain and filtering saturated with energy, and loss saturated with power. It is found that 
general initial conditions evolve into dark solitons under appropriate requirements also met in the experimental 
observations. The resulting pulses are well approximated by dark solitons of the unperturbed nonlinear 
Schrodinger equation. Notably, the same framework also describes bright pulses in anomalous and normally 
dispersive lasers. © 2010 Optical Society of America 
OCIS codes: 140.4050, 140.3510, 060.5530, 190.5530 



Laser frequency stabilization via mode-locked (ML) 
lasers has become an indispensable tool in many re- 
search activities. Advances in optical frequency stan- 
dards have resulted in the development of precise fre- 
quency measurement capability in the visible and near- 
infrared spectral regions. Although the potential for us- 
ing ML lasers in optical frequency synthesis was rec- 
ognized early, the available lasers did not provide the 
properties necessary for fulfilling this potential until re- 
cently [T]. The recent explosion of measurements based 
on ML lasers has been largely due to the development of 
the Kerr-lens mode-locked ti:sapphire laser and its capa- 
bility to generate sufficiently short pulses. Alternatively, 
amplified pulses of a similar duration are created through 
the compression of pulses via self-phase modulation in a 
hollow core fiber [5]. 

Dark solitons are intensity dips on a constant back- 
ground with a phase jump across their intensity min- 
imum. Since their discovery, cf. j^, dark solitons have 
attracted considerable attention, especially in the fields 
of nonlinear optics j4j and Bose-Einstein condensates [5] . 
In the non-mode-locking regime, the first train of dark 
solitons was successfully achieved in an all normal dis- 
persion erbium-doped fiber laser with an in-cavity po- 
larizer. In ML lasers they have been difficult to gener- 
ate with the first experimental observation only recently 
reported [2]; see also [S]. Apart from the bright pulse 
emission, it was observed that, under under strong con- 
tinuous wave emission, appropriate pump strength and 
negative cavity feedback, a fiber laser can also emit single 
or multiple dark pulses [2]. 

To analyze dark solitons in ML lasers we use a 
model based on a variant of the nonlinear Schrodinger 
(NLS) equation, termed the power-energy saturation 
(PES) equation [7], suitably generalized to include non- 
vanishing boundary conditions at infinity. It is expressed 



in the following dimensionless form. 
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where the complex electric field envelope ip{z,t) is sub- 
ject to the boundary conditions \'il}{z,t)\ — > iV'ool as 
\t\ oo. Here, E{z) = /^^(iV'ooP - IV'I^) dt is the dark- 
pulse energy, P{z,t) = IV'ooP — IV'P is the instantaneous 
power, while Eq and Pq are related to the saturation 
energy and power, respectively. Furthermore, g, r, and 
/ are all positive, real constants, with the correspond- 
ing terms representing saturable gain, spectral filtering, 
and saturable loss. Typical dimensional numbers can be 
found in Ref. 0. 

The above definition of the pulse energy, although un- 
conventional, corresponds to the physical properties of 
the system. Indeed, as the peak-to-background intensity 
ratio can achieve large values in bright pulses, so can 
the focused-to-unfocused vacuum ratio in dark pulses. 
As such, dark pulses can be thought of as focusing the 
vacuum [9] . The definition of power follows consistently 
from E = P dt. A modified model (yielding similar 
results), with more standard definitions of energy and 
power, is also discussed below. 

To study the dynamics and mode-locking capabili- 
ties of our model, we integrate Eq. ([T]) with a fourth- 
order Runge-Kutta method, with initial profile "0(0, t) = 
tanht, though others, with a 7r-phase jump, can be 
used as well. The gain parameter g is varied, while 
Eq = Pq = 1 and t = Z = 0.1. In the top panel of Fig. 
[1] we show the evolution of the background amplitude, 
which sets the dark soliton amplitude, for different values 
of the gain parameter g. Locking onto stable dark soli- 
tons is only achieved when the gain term is sufficiently 
strong, i.e. the parameter g is large enough to counter 
balance the losses. Thus, for g ^ OA (see dashed line in 
the top panel of Fig. [T]), the dark soliton decays quickly 
due to excessive loss. On the other hand, for g = 0.2, the 
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Fig. 1. (Color online) Top panel: evolution of the pulse 
background of an arbitrary initial profile under the PES 
equation with different values of gain. Bottom panel: 
complete evolution of the initial profile for g ~ 0.5. 

pulse initially decreases, but then locks to a specific am- 
plitude and width. Finally, for g = 0.5, 0.7, 1.0, a stable 
evolution characterized by an increase of the soliton am- 
plitude is obtained almost instantaneously; see bottom 
panel of Fig. [T] 

To determine the evolution of the background wave 
in the framework of Eq. ([l}, we assume that ipi^) = 
ipooiz) expli6{z)]. Separating real and imaginary parts, 
yields the following equation for the background ampli- 
tude iJoo{z): 



dV'o 
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Here, an approximate solution of Eq. ([l} in the form 
\^p(z,t)\ = l^/itxj tanh{ipoot)\ is assumed, which gives E = 
2 1 V'oo I - Equation Q also describes the evolution of V'oo, 
as depicted in the top panel of Fig. [l] and illustrates the 
importance of energy saturation. Indeed, if we consider 
the so-called similariton supporting equation (NLS with 
linear gain) corresponding to t — I = [SlllOj. with en- 
ergy saturation absent (2|'(/;oo|/£'0 = 0), the pulse ampli- 
tude grows exponentially at a rate defined by the gain g. 
Contrary, when saturation is present, this exponentially 
growing state cannot be reached. If fact, all terms in Eq. 
p]) are important for successful mode-locking. Spectral 
filtering is important in maintaining a short pulse dura- 
tion, which is a generic phenomenon to mode-locking in 



the normal dispersion regime [TT]. With spectral filter- 
ing in the cavity, stable ultrashort pulses can be gener- 
ated at very high values of the group- velocity dispersion. 
Saturable loss results in higher net gain for a train of 
short pulses (compared to continuous- wave operation), 
which is also necessary for mode-locking. In some mod- 
els, loss is introduced in the form of fast saturable power 
absorbers placed periodically [12]. It has been shown 
that the lumped and distributive models yield similar 
results [TB] . 

A stable equilibrium (attractor) exists for Eq. ([2]) and 
can be found setting di/'oo/d^ = 0, namely. 
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Equation 1^ is the resulting background amplitude of 
the dark soliton and agrees with direct numerical simu- 
lation, as seen in Fig. (3] Thus, dark solitons tend to an 
equilibrium (mode-lock) with constant energy and back- 
ground. 

The above findings agree with the experiment of 
Ref . [2] , where the mode- locked dark pulses were identi- 
fied as NLS dark solitons. In particular, they are black 
solitons (stationary kinks) characterized by a 7r-phase 
jump across the soliton notch and they can be de- 
scribed analytically by a hyperbolic tangent profile, i.e., 
ipiz, t) = V'oo tanh('!/'ooi)- In Fig. [5] we plot the black soli- 
tons of the PES and NLS equations (for the same value 
of V'oo, as determined from Eq. ([3])) for different values 
of g. The amplitudes match so closely that they are in- 
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Fig. 2. (Color online) Black solitons of the PES and NLS 
equations for different values of g. 

distinguishable in the figure, meaning that the combined 
perturbation of gain, filtering and loss in the PES model 
underlies the mode-locking mechanism: its effect is to 
mode-lock to a black soliton of the pure NLS equation 
with the appropriate background, Eq. (j3|). On the other 
hand, mode-locking to grey solitons, namely moving dark 
pulses with a phase-jump less than tt at their (non-zero) 
intensity minima, are not found in this model as they 
have insufficient dark-pulse energy. This too agrees with 
Ref. [2] where grey solitons were not observed. Note that 
recent experimental [M] and theoretical [15] results in 
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ML lasers in the anomalous dispersion regime indicate 
that the normalized intensity of a pulse in a ML laser can 
be described by the bright solitons of the unperturbed 
NLS equation. 

The conventional definition of the pulse energy may 
also be used in Eq. ([T]), namely E = j: /Tt'/2 IV'P dt, 
where T is the averaged time through the cavity. In this 
case, looking for stationary solutions of Eq. @ yields 
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for the background where, for sufhciently large T, E ~ 
T /-T/2l^l^ « i|VoopT = IV-ooP- Equation (gD is 
a more restrictive condition for mode-locking which can 
also be attributed to the experimental difficulties of the 
problem. Furthermore, this definition results in more 
complicated dynamics and the development of a shelf 
on the soliton background, which is more pronounced as 
compared to the previous case. In fact, as shown else- 
where |16| . shelves occur naturally in the perturbation 
theory for dark solitons. We illustrate the resulting shelf 
in Fig.|3l where now g ~ 0.2, 1 ~ 0.1, Eq = 1 and Pq = 5. 
Due to the noisy background this feature maybe difficult 
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Fig. 3. (Color online) The development of a shelf in the 
PES equation (solid line) as compared to the solution of 
the NLS equation (dashed line). 

to observe in an experiment. However, the shelf also af- 
fects the phase of the resulting pulse (see inset in Fig. [3]) 
and since the experimentally observed pulses exhibit a 
sharp TT-phase jump |2], the first dark-pulse energy defi- 
nition may give a closer approximation for this problem. 

Lastly, we mention the notion of multiple dark soli- 
tons. To observe these, we evolve Eq. ([1]) with initial 
condition ■0(0, i) = V'oo tanh('!/)oo^ + ^o) for t < and 
■0(0, t) = —ipoo tanh(-0tx)i — ^o) for t > 0, where, tpoo can 
be approximated by ipoc/'^] the dark-pulse energy is now 
twice that of the single soliton. In Fig. S] we show its 
evolution. As in the case of bright pulses [T7], the two 
solitons are taken initially sufficiently far apart (about 
10 times their full- width of half-maximum) to minimize 
interactions. The major difference with the bright case, 
where the interaction is logarithmically slow [17] , is that 
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Fig. 4. (Color online) Two-dark-soliton evolution. 



pulses repel sooner here due to the shelf interactions. 
While the solitons are stable for reasonably long distance 
enough to also be experimentally observed, nevertheless 
one can expect this case to be sensitive as with the single 
dark soliton. 

In conclusion, we presented a power-energy saturation 
model that describes mode-locking of dark solitons in 
lasers. Much like their bright counterparts, these pulses 
(as well as multiple pulses) exist when sufficient gain 
is present in the system, which agrees with the experi- 
mental observations. The specific energy saturated gain 
and filtering, and power saturated loss, are crucial to the 
mode-locking mechanism. The resulting pulses are essen- 
tially modes of the unperturbed NLS equation with back- 
ground amplitudes appropriately defined by the gain and 
loss parameters. 
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Force Office of Scientific Research, under grant FA9550- 
09-1-0250. 
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